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ABSTRACT 

The Gamma-Ray burst (GRB) - afterglow transition is one of the most interesting and 
least studied GRB phases. During this phase the relativistic ejecta begins interacting 
£T) . with the surrounding matter. A strong short lived reverse shock propagates into the 

' ejecta (provided that it is baryonic) while the forward shock begins to shape the sur- 

rounding matter into a Blandford-McKee profile. We suggest a parametrization of the 
early afterglow light curve and we calculate (analytically and numerically) the observed 
parameters that results from a reverse shock emission (in an interstellar medium [ISM] 
environment). We present a new fingerprint of the reverse shock emission that is added 
Tj" , to the well known t~ 2 optical decay. Observation of this signature would indicate that 

■ the reverse shock dominates the emission during the early afterglow. The existence 

. of a reverse shock will in turn imply that the relativistic ejecta contains a significant 

(*S baryonic component. This signature would also imply that the surrounding medium 

is an ISM. We further show that: (i) The reverse shock optical flash depends strongly 
' on initial conditions of the relativistic ejecta. (ii) Previous calculations have generally 

\ overestimated the strength of this optical flash, (iii) If the reverse shock dominates 

the optical flash then detailed observations of the early afterglow light curve would 
possibly enable us to determine the initial physical conditions within the relativistic 
ejecta and specifically to estimate its Lorentz factor and its width. 

Key words: gamma-rays: bursts-shock waves-hydrodynamics 



1 INTRODUCTION 

According to the internal-external shocks model (Piran & Sari, 1998) the prompt gamma-ray burst (GRB) is produced 
by internal shocks within a relativistic flow while the afterglow is produced by external shocks between this flow and the 
surrounding matter. The early afterglow appears during the transition from the prompt 7-ray emission to the afterglow. 
During this transition the relativistic flow, ejected by the source, interacts directly with the circum burst medium. This 
interaction can be used to pin down the nature of the relativistic flow (baryonic or Poynting flux). In a baryonic flow the 
reverse shock (RS) that propagates into the ejecta produces both optical and radio emission. With Poynting flux we expect 
only the higher energy forward shock (FS) emission. While other sources of early optical and radio emission may exist also in 
Poynting Flux flow, we show here that the RS emission has a very robust optical and radio observable signatures that is very 
unlikely to be imitated by other phenomena. If the flow is found to be baryonic then the early afterglow signal could serve as 
a diagnostic tool for the properties of the ejecta. This in turn, would shed light on the nature of the inner engine that powers 
the GRB. 

Numerous authors (Meszaros & Rees 1997; Sari & Piran 1999a; Sari & Meszaros 2000; Soderberg & Ramirez-Ruiz 
2003a; Zhang, Kobayashi & Meszaros 2003) considered the emission from the reverse shock. Strong optical flashes in a rough 
agreement with the RS predictions (Sari & Piran 1999b; Meszaros & Rees 1999; Wang, Dai & Lu 2000; Fan et al. 2002; 
Soderberg & Ramirez-Ruiz 2003b; Fox et al. 2003a; Weidong et al. 2003; Kumar & Panaitescu 2003) were observed in two 
bursts (GRBs 990123 and 021211). In the first the RS predicted radio flare was observed as well (Sari & Piran 1999b; Kulkarni 
et al. 1999). On the other hand the early (the first lOminutes) optical emission observed in two other bursts (GRBs 021004 
and 030418; Fox et al . 2003b; Rykoff et al. 2004) did not agree with the simple predictions of an RS emission (see however 
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Kobayashi & Zhang 2003a). Furthermore, upper limits of ~ 15th mag on the prompt optical flux of several bursts (Williams 
et al. 1999; Rykoff & Smith, 2002; Klotz et al., 2003; Torii 2003a,b) have lead to the so called "optical flash problem": lack 
of bright optical flashes (corresponding to RS emission) in many bursts. 

Swift is expected to provide a large number of deep (~20 mag) early (~minute) optical observations. We provide here 
detailed predictions of the RS emission of a baryonic flow that interacts with a constant density circum-burst medium (such 
as an interstellar medium [ISM]). These predictions can be confronted with the upcoming observations. We show that when 
the early afterglow is found to be dominated by such RS emission (by passing the observational tests) then the observations 
enable us to determine the initial physical properties of the relativistic outflow and to constrain the microscopic parameters 
in the emitting region. We also show that the peak flux depends sensitively on the strength of the reverse shock. It can 
vary over more than five optical magnitudes between a mildly relativistic and ultra relativistic shocks. Furthermore, previous 
calculations have generally overestimated the peak flux of the mildly relativistic RS by up to 7 optical magnitudes! In fact 
the FS emission may even dominate over the RS emission at early times. Altogether these results suggest a solution to the 
"optical flash problem" . The calculated radio light curve shows that the radio flare lasts long after the optical flash. We find 
that in the generic case, in addition to the typical decay of the optical flash (Sari & Piran 1999a), the flash to flare time ratio 
and intensity ratio provide another new test that the emission results from a reverse shock. 

Other mechanisms apart from the RS in an ISM environment can produce an optical flash. Few examples are internal 
shocks (Meszaros & Rees 1997), a pair avalanche (Thompson & Madau 2000, Beloborodov 2002), RS produced in a wind 
environment (Chevalier & Li 2000; Wu et al., 2003; Kobayashi & Zhang 2003; Kobayashi, Meszaros & Zhang 2004) and 
even the FS as we discuss here. While these mechanisms might be able to produce a bright early optical flash, they are not 
expected to produce the combination of the early optical and radio emission that we show here to arise from an RS in an 
ISM environment . Thus observing this signature in many bursts in the future will solve further open questions than just 
the flow's constitution (Baryonic vs. Poynting flux). It will also reveal the circum-burst density profile and determine the 
dominant mechanism that contributes to the optical flash. In a separate paper (Nakar & Piran, 04) we apply the tests and 
the diagnostic tools presented here to GRB 990123. We find that its early afterglow emission is remarkably consistent with 
our predictions of the RS emission, suggesting strongly that at least is this case the early afterglow is produced by a baryonic 
flow propagating into an ISM. 

The structure of the paper is somewhat unusual. The RS emission in the most general case presents a large and complex 
diversity, but the generic behavior is rather simple. Therefore, we begin in 33 with a summary of the generic optical and radio 
observables, and demonstrate how to use them in order to confirm that the emission results from an RS and in this case 
to determine the physical properties of the relativistic outflow. Later we describe the general analytical theory (iJSJ and the 
numerical simulations (2J. 



2 THE GENERIC EARLY AFTERGLOW LIGHT CURVE 

Consider a homogenous 2 cold baryonic shell expanding relativistically into an homogenous cold inter-stellar medium (ISM). 
The problem is well defined by the shell's (isotropic equivalent) energy E, width A, initial Lorentz factor r o and the ISM 
density n. As the ejecta shovels the ISM, a forward shock and a reverse shock are produced. The nature of the RS is determined 
by the dimensionless parameter £ = (Z/A) 1/2 r~ 4/3 (Sari & Piran 1995, hereafter SP95), where I = (3E/(ATmm p c 2 )) 1/3 . If 
£ <C 1 the RS is relativistic and most of the shell's bulk motion energy is dissipated in a single shock crossing of the shell, 
which occur at a radius -Ra ~ Z^A 1 / 4 (SP95). For £ S> 1 the RS is Newtonian and many crossings are required to dissipate a 
significant fraction of the energy. In the generic case £ < 1. In this case after the RS crosses the ejecta once, the circum-burst 
gas, shocked by the FS, forms a Blandford-Mckee profile and the original ejecta expands and cools down at the tail of this 
profile (Kobayashi & Sari 2000; hereafter KS00). Therefore the hydrodynamical (and the radiative) evolution is separated to 
two phases - during the RS and after the RS. The emission from the separation point reaches the observer simultaneously 
with a photon emitted on the line of sight from the rear end of the shell at Ra'- 



where c is the light speed, z is the redshift and Mt = 1.4. Here and through out the paper A4 stands for numerical correction 
factors to the analytic estimates (see secll). 

The evolution before to is highly sensitive to the initial profile of the shell and in particular to the value of f . Thus the 
observables before and at to depend strongly on the initial properties of the shell, and as such can be used as a diagnostic 
tool of these properties. On the other hand, the RS that crosses the shell erases, to a large extend, the initial shell profile. 

1 For example, in the pair avalanche process the optical flash results from a pair enriched FS and this emission is becoming harder with 
time (Beloborodov 2002), in contrast to the RS emission. Thus, no correlated radio flare is expected. 

2 The limits of the results presented in this section when the shell is inhomogeneous are discussed in sec. [3] 
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Moreover, the evolution during the expanding and cooling phase depends only weakly on £ (KSOO). Therefore, the behavior 
after to is insensitive to the initial conditions and as such it provides a very unique and easily identified signature of an RS 
emission. 

The general behavior of a basic RS (i.e. with no complications as refreshed shocks during the RS etc.) over a broad 
spectrum is described in Tabled In this section we summarize only the generic behavior of the optical and radio observables 
and demonstrate how to determine the physical parameters, £, A, and F D , from the optical light curve. This should be used 
with care as a non-generic behavior is always a possibility. 

The RS optical flash peaks at to. Therefore, similarity to late afterglow parametrization (Beuermann et al. 1999), we 
suggest to parameterize the RS optical emission as: 

where cti > and «2 < are the power-laws indices of a broken power-law that peaks at 3 ~ F£ opt (to) = Fq. s determines 
the sharpness of the peak. The optical frequency, v ovt is expected to satisfy v^^m < v o V t < v r c where v r a ,v r m and v r c are the 
self-absorbtion, synchrotron and cooling frequencies in the RS respectively. In this case (Sari & Piran, 1999a): 

a 2 « -2. (3) 

The decay slope, 012, as a post to observable, is very robust. KSOO show numerically that oli w —2 for various values of £ and 
p and thus it is a signature of a generic RS emission. On the other hand, ax is most sensitive to £. When £ <C 1, ai ~ 0.5 and 
as £ increases so does ot\. For 0.05 < £ < 5 it can be well approximated as (see Fig. Ep.) : 

cti ~A/; :1 (0.5 + !(e-0.07£ 2 )), (4) 

where p is the power-law index of the electrons' energy distribution, and Af a ,i = 1.2 . Thus, a measurement of a\ can determine 
the value of £ (up to the uncertainty in the value of p). Once £ is known, Eq. Q, is solved for A. Having £ and A one can 
find r o using: 

r o = 188r 3/4 Ar 2 3/8 (iWn) 1/8 , (5) 

where we denote by Q x the value of the quantity Q in units of 10 x (c.g.s). Note that when ^ < 1 both Eqs. Q and 
are insensitive to £ and only a lower limit of F can be found. r o depends very weakly on the ratio E/n. Finally we find 
numerically that the sharpness parameter s depends strongly on the initial profile of the shell, but not on £. The larger the 
initial Lorentz factor dispersion (5T = r o , m0 i/r o , m m) is the smaller is s (wider peak). A homogenous shell (8T = 1) results 
in a very sharp peak, s ~ 10, while mild dispersion of ST — 2 may be sufficient to reduce s to w 1. 

It is remarkable that these initial parameters can be determined without using Fq, and thus with no dependance on the 
poorly known internal parameters, e e and es. The value of Fq can be used to constrain these parameters: 

-±b£, »2B_« /sfo-aA 1 ' - ' 1 



FS = 0.6mJy A/>(1 + z)-*? 1.5 2 5 " p (^i)' 



P+l „_i_2 i _|_ P 3p 



(6) 



where the numerical correction factor is Mf ~ 1/5 and all the parameters and notations are as in Tabled The function 
A r FQ {^) is approximated in the range of 0.1 < £ < 2.5 by: 

A r F ,o(0 « 180£°- e5 (6 ■ 10~ 4 r 2 - 6 ) ^ • (7) 

The exact value of A r F is given in Eq. 1131 . and must be used when £ is outside of the range above. Fg depends strongly 
on £ and it varies by 2 orders of magnitude within the range most relevant for GRBs (0.1 < £ < 3). The relativistic and the 
Newtonian approximations (Fig.^j) overestimate Fg. Specifically the commonly used Newtonian approximation overestimate 
F£(£ = 1) by a factor of 200. The numerical correction factor, Mf, adds another factor of 5! 

The radio emission continues to rise after to and it peaks at a later time, t,, when v ra dio = v r a . This happens during 
the "cooling" phase of the shocked shell material. Therefore the radio behavior both before and after t» is a robust feature 
(i.e. does not depend on the initial conditions). The light curve depends only on the relations between v T m , v r a and v ra dio and 
the only remaining influence of the initial conditions is via the values of the break frequencies at to- Over a wide range of £ 
values i4(t ) < i£(to) ~ 10 12_13 Hz 4 . In this case the radio flux at t > to can be also characterized by the parametrization of 



3 For s > 1 the peak of Eq. J5J is obtained at w to and its value is w Fq. Numerical simulations show that in the case analyzed here 
indeed s J> 1 (see sec. Thus, to and Fq that are found according to the best fit of Eq. can be taken, for any practical purpose, 
directly from the observations as the time and the Flux of optical peak. 

4 A different case than the generic one(i.e different relation between v^, i/^ and is r adio)is more likely here than in the optical emission, 
specially when the RS is ultra-relativistic. 
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Eq. © with (see Table [TJ 5 : 



a,.,i » 


1.25 


t. 




to 




F, 




T 




Eq. JHJ predicts a 


that the 


emission 



Qr,2 ~ —2 
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F* fu\ 2 + 3 _ / v opt 



,, i , , . !000. (9) 

Note that this value can be larger or smaller by a factor of ~ 3 (for a given p), because of the uncertainty in the post-RS 
dynamics (KSOO). Together with the optical decay signature, Eqs. I8I9H provide a unique imprint of a baryonic RS. The 
determination of fa(to) (and maybe even i/^(to)) provide additional constraints on e e and eb. 

The early afterglow behavior described above is relevant only when the RS is produced by interaction with an ISM like 
densities or lower (n < 100cm -3 )). This condition is not satisfied if the circum-burst medium is the wind of a massive star. 
For any reasonable parameters of such a wind the external density during the crossing of the RS is few orders of magnitude 
larger than this of a typical ISM. This brings the cooling frequency well below the optical bands (Chevalier & Li 2000) and 
possibly the self absorbtion frequency above the optical band (Kobayashi et al. 2004) . This changes of the frequencies sequence 
changes also the resulting behavior of the optical light curve (e.g. the value of 02) . Therefore observing an early afterglow 
that shows the RS emission signature described above implies also that the density of the external medium is an ISM like. 



3 THEORY 

The nature of the RS is determined by the dimensionless parameter £, which in turn determines the ratio a, of the Lorentz 
factor of the shocked matter (in the explosion rest frame), y T , to r o : 

a = lr /T . (10) 

a can be derived directly from the relativistic jump conditions (SP95). It satisfies 6 : 

(12/£ 3 - l)a 4 + 0.5a 3 + a 2 + 0.5a -1=0. (11) 

In the relativistic regime a ~ an = 12 _1 ' 4 £ 3 ' 4 , while in the Newtonian regime a » ajv = 1. Both approximations overestimate 
a in the intermediate regime and the deficiency is largest when £ « 1 (Fig. JQi)). For f < a few the RS emission peaks when 
the RS reaches the back of the shell. At this stage the pressure, p r , and the density, n r , in the shocked shell as measured in 
the shocked fluid rest frame, are: 

Pr = ^a 2 V 2 nm p c 2 ; n r = fnT 2 (2(a + l/o)/3 + 1). (12) 

Assuming a homogeneous initial shell and homogenous conditions within the shocked region 7 these hydrodynamical conditions 
determine v r a ,^ m , v r c and the peak flux F^ max at t — to- These values appear in Tableland can be used to estimate the 
radio and the optical emission at to- Using these values we derive the optical emission, Eq. @ and the exact value of A r F (£) : 

A r F , (0 = 770a 4p-2 (a 2 + 1.5a + i)1-p^-3.75 P(1 + 0.5^!)^. (13) 

Next we consider the evolution at t < to- The flux at t < to can be determined by parameterizing all the quantities 
according to the fraction, /, of the shell that the RS has crossed: A(/) and E(f) on f while n is constant. This implies 
£(/) « /~ 1/3 and R(f) tx / 1/2 . The observer time 

t(/)cx/(l + 0.5Me(/) 3/2 ) (14) 

and the optical flux (for v r a , v m < v ovt < v T c ; see Table 0: 

F r v , 01 M) « fa 4p - 2 (a 2 + 1.5a + l) 1 "*, (15) 

combine to yield ai. In the relativistic regime (a tx £ 3//4 oc / _1,/4 ) t oc / and Fo oc f 1 ' 2 , hence a\ = 0.5. When £ increase 

5 As long as v ra dio < "mi a r ~ 0.5. Observing the transition time to a r> i fs 1.25 determines i/J^to). The early behavior at t < to can 
be found from Table HI 

6 For a relativistic adiabatic constant - 4/3. Our results do not change significantly if the adiabatic constant varies smoothly between 
4/3 when f < 1 to 5/3 when f > 1 

7 These assumptions are relaxed in the numerical simulations 
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Figure 1. The parameter, a and various break frequencies and flux densities as functions of £. (a): The parameter a according to Eq. 
Hill (thick line), and its relativistic (dashed line) and Newtonian (thin line) approximations, (b): Fq according to Eq. IS1 fc 1131 when 
a is calculated according to Eq. 1111 (thick line), and when the relativistic (dashed line) or Newtonian (thin line) approximation of a is 
taken. The physical parameters considered here are similar to these considered in panel d. (c): The ratio between the RS peak flux and 
the FS flux at the same time. e e and eg are assumed to be similar in the RS and the FS. (d): The synchrotron (solid line) and the self 
absorbtion (dashed line) frequencies of the RS emission at to- The parameters considered here are: E52 = 1, e e = 0.1, eg = 0.01, n = 1, 
p = 2.5, t = lOOsec and 2 = 1 (D 2S = 1). 



the logarithmic time derivative (dlog(F^ opt ) / d\og(t)) varies with time, and its value at t < to depends strongly on £. The 
description of the light curve as a power-law with an index ct\ is only an approximation. We estimate cvi as the mean value of 
dlog(F^ opt )/dlog(t) during to/2 < t < to, and compare it to the standard deviation in this time range (Fig. The small 
deviation compared to the mean value justifies the power-law approximation. 

The evolution at t > to is dictated by the post-RS hydrodynamics. This hydrodynamics evolution was investigated both 
analytically and numerically in KS00, and found to be almost independent of £. We use the hydrodynamic evolution presented 
in KS00 to determine the optical and radio light curves. 

Finally we calculate the contribution of the forward shock emission. The FS hydrodynamical conditions at Ra are: 

If = It ; Pf = p r ; rif = AaY n. (16) 

The corresponding spectral parameters at to, {vl,m,c(pa) an d F^ max (to)), are listed in Table □ The FS emission, in contrast 
to the RS emission, depends rather weakly on £. Thus, the ratio between the two vary strongly with £. Specifically, the widely 
used relation v^/v^ = (n r /nf) 2 depends strongly on £. For £ « 1 it indeed equals » r„. However, for £ = 0.1 this ratio equals 
~ 5 ■ 10~ 4 r„. Similarly, the ratio F^ opt (to) / F^ opt (to) determines whether the RS or the FS dominates the emission. Fig. 0; 
depicts this ratio for typical parameters. With these specific parameters, and assuming similar e e and cb in both regions, the 
FS optical emission cannot be neglected for £ > 0.5. 



4 NUMERICAL SIMULATIONS AND INHOMOGENEOUS SHELLS 

The analytic calculations presented above include several approximations. In order to verify the accuracy of the analytic 
calculations we carried out detailed numerical simulations (Nakar & Granot 2004) of the early afterglow emission. We use 
these simulations to determine numerical correction factors, denoted Mx, for the analytic calculations. The hydrodynamics 
simulations were done using a one dimensional Lagrangian code that was provided to us generously by Re'em Sari and Shiho 
Kobayashi (KS00). The synchrotron radiation code is described in Nakar & Granot (2004). This code provides an accurate 
synchrotron light-curve and spectrum, taking into account the realistic hydrodynamical profile 8 of the emitting region, the 
exact heating and cooling history of the electrons and the precise photons arrival time to the observer from each radiating 
clement. 

We have carried the simulations for a range of parameters with 0.05 ^ £ ^ 3. Using these simulations we obtain 
numerical corrections coefficients and determine the accuracy of the analytic estimates of to, Fq , F^ opt and a\. We find that 
when including the numerical correction factors (that range from 0.2 to 1.4), Eq.Q for to is accurate up to 10% while Eq.© 
(for F^ opt ) and the expression for F^ opt are accurate up to a factor 2. The sharpness parameter s was not considered before, 
and we find its value numerically. We find that if the shell is homogenous the peak is very sharp, s ~ 10, regardless of the 

8 We neglect the feedback of the radiation energy loses on the hydrodynamics. This is justified in the likely case that t e «l 
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Figure 2. Results of detailed hydro+synchrotron numerical simulations, (a): a\ (£) as estimated by the mean value of dlog(F^ t ) /dlog(t) 
al to/2 < / < to, according to Eqs 1111151 (thin line), according to Eq. Q) (dashed line) and according to the numerical simulations 
(dots). The standard deviation of dlog(F£ opt )/dlog(t) at to/2 < t < to is depicted by the shaded area (Eqs 1141151 and the error bars 
(numerical simulations), (b): Five numerical (normalized) optical light curves with various £ values, and their best fits according to Eq 
(see sec. [I). 



value of £. Fig. [5Jj depicts five numerical light curves with 0.05 ^ £ ^ 2, and their best fits according to Eq.0. In all these 
fits Fq is within factor of 2 of Eq.@, to is within 10% of Eq.Q, ot\ is within the spread of Eq.@ (the shaded area in Fig. 
Hi), a 2 « -2 and s « 10. 

So far we have considered homogeneous shells. Clearly, the light curve resulting from an inhomogeneous shell would 
depend on the shell's profile. To investigate partially the effect of inhomogeneity, we have carried out numerical simulations 
of shells with linear Lorentz factor profile (ST = 1.6) and a constant energy per a rest frame length interval. The value of £ 
varied between 0.05 and 1 where F is taken as the mean value of the initial Lorentz factor. As expected in all the cases we find 
Q2 ~ —2 (which indicate that also the radio emission should be insensitive to the shell's initial profile) and a sharp decrease 
in s when £ » 1 (s « 1 compared to s w 10 in the homogenous case). Fq , to and Qi are similar to the values obtained in the 
homogenous case (a maximal difference of a factor of 1.5 between the homogeneous and inhomogeneous cases) 9 . These results 
make us confidant that the solution we present here for homogenous shell is generic and applicable also for inhomogeneous 
shells (at least as long as 8T is not much greater than 1, and the Lorentz factor profile rises monotonically and regularly). 

Finally we discuss briefly two caveats, which may arise naturally in a baryonic relativistic flow, and may alter the RS 
emission and produce a non-generic RS optical light curve. The first caveat is a slow tail of the wind with a small but not 
negligible, amount of energy (compared to the total energy of the wind). Such a tail may results from the adiabatic cooling 
of the wind after the phase of the prompt 7-rays emission, during which the wind must be hot. In this case the peak may 
be obtained before the RS finishes crossing the shell. Therefore the decay after the peak is expected to be shallower then 
the generic value of —2 (even for several orders of magnitude in time), becoming gradually steeper until the generic value is 
obtained 10 . The second caveat is a highly irregular density profile. Such profile may result from a highly irregular ejection 
process in the source, as expected in the internal shocks scenario when a burst is highly variable. Since a hydrodynamic 
evolution smoothes the pressure and the velocity profiles, but not the density profile, such irregularities are expected to be 
carried by the flow also to the RS phase. If these irregularities in the density profile are large enough, they are expected to 
be reflected in the RS optical light curve during its rising phase and before its decay reaches the asymptotical value. Thus, 
a detailed observations of this phase may reveal the exact profile of the ejecta, and maybe even be used to test the internal 
shocks model (as the density profile in the flow is expected to reflect the light curve produced by internal shocks). Practically, 
if the early afterglow light curve is highly irregular, with no underlying power-law, then the analysis method described here is 
not applicable, and a theory that describes highly inhomogeneous RS is required. Note, however, that if the asymptotic value 
of the decay is reached soon after the peak, the tests of the RS emission (Eqs. |21&|5J are still applicable. 



9 Note that with a small value for s (broad peak) qi reaches its "asymptotic" value only far from the peak. 

10 This may explain the shallow early decay of GRB 021004. This scenario is very similar to the continuous refreshed RS introduced by 
Sari & Mcszaros 2000, and suggested by Fox et al. 2003b to explain GRB 021004. The only difference is that here the slow tail of the 
flow is produced naturally from the hydrodynamics of the relativistic wind and not by a special source activity. 
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Table 1. The Break Frequencies and Maximal Flux Densities 
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Notations - e e , £3: the fraction of the internal energy in relativistic electrons and magnetic field respectively; p: the electrons spectral 

index; D: proper distance to the burst; e e = 30e e (p - 2)/(p - 1); £> c = (1 + /V t 0.5£ 3 / 2 ); ha S a 2 + 1.5a + 1; = A£(f)M£(l) 

Using the table - F u (v) is found according the maximal flux, F v ,max, the values of the break frequencies and the spectral power law 

indices between them. All these vary between the different cases, where at each case F V} max is at a different break frequency: 

v m <v a < v c (generic case): F^ t , nax (u = u a ) = F Vtb (v a / V m )( 1_p )/ 2 : ; Power law indices: 2 < v m < 2.5 < v a < (1 - p)/2 < v c < —p/2 

v a < Vm < v c : F Vtmax (y = Vm) = F^b ; Power law indices: 2 < v a < 1/3 < v m < (1 - p)/2 < v c < -p/2 

v a <v c < v m '- F Vt max(v = v c ) = F v<h ; Power law indices: 2 < u a < 1/3 < v c < —1/2 < v m < — p/2 

Column (2): The values at t = to- Column (3): The dependance at t ^ to on £. Column (4): The evolution at f < to can be found 
using this column and Eq. 1141 (see sec.|2J. Column (5): The approximated power-law indices at t > to, these values are uncertain by 
at least ±0.05 — 0.1 (see KS00). The values are calculated using the hydrodynamics of Eqs. {JUEU (t < t ) and KS00 (t > t ), and 
the radiation calculations described in Sari, Piran & Narayan (1998), Granot, Piran & Sari (1999) and Granot & Sari (2001). 



*v[ = v r c at f < to 



Va < V m < V c 



' V a < V c < Vn 



ttt 



Vm < V a < V c 



5 CONCLUSIONS 

We have suggested here a parametrization of the early optical emission, as a broken power law with five parameters. We have 
calculated the values of these parameters for a RS produced by the interaction of baryonic wind with a circum burst ISM. 
Our main conclusions are: (i) The optical decay, {0.2 ~ —2) and the consistency between the peak time and flux of the optical 
and the radio (Eq. |HJ are robust features of a RS emission in an ISM environment over a large range of initial parameters. 
Observations of an early optical emission with these features would suggest that: a) A significant fraction of outflow energy is 
baryonic 11 . b) The circum-burst medium is an ISM like, c) The RS emission is dominant over other possible sources of optical 
flash, (ii) The values of the observables before the optical peak depend strongly on the strength of the shock, and can be 
used to pin down the initial conditions of the flow, (iii) The combination of optical flash and radio flare may also constrain 
the microscopic parameters in the emitting region. 

In addition to the specific optical predictions we presented detailed analytical results of the expected emission over the 
whole spectrum. The advantages of these calculation over previous ones are that they do not make any approximation on the 
strength of the RS (i.e. relativistic or Newtonian) , and that they are confirmed (and corrected) by numerical simulations. The 
conclusions of these calculations are: (i) Previous calculations overestimated the intensity of the optical flash. Most pronounced 
is the Newtonian approximation for £ w 1 that overestimate the optical flash by up to three orders of magnitude, (ii) An 
optical flash brighter than 15th mag is expected in some but not in all GRBs. A GRB with typical parameters and moderate 
energy = 1), is expected to produce a maximal flux of R m ag ~ 17 — 19 when the RS is mildly relativistic (£ « 1). (iii) 
Over some reasonable range of the parameters space the FS emission dominates at all times, (iv) When the RS is relativistic 
(£ <C 1) most of the emission is released in the optical. When £ increase, the emission is shifted to lower energy bands but it 
does not reaches the radio, as v a in this case is ~ 10 12-13 Hz. These results suggest a solution to the "lack" of optical flashes. 
Only a fraction of the flashes is expected to be bright enough for detection in the current observations, and in some cases an 
FS emission is expected to dominate from the beginning. 

We thank R. Sari and S. Kobayashi for providing us with their relativistic hydrodynamics code and J. Granot, P. Kumar, 
R. Mochkovitch, F. Daigne and E. Rossi for helpful discussions. The research was supported by the US-Israel BSF and by 
EU-RTN: GRBs - Enigma and a Tool. EN is supported by the Horowitz foundation and by a Dan David Prize Scholarship 
2003. 



Clearly a small fraction of the outflow energy as a Poynting flux would not change the results presented here, while a very large 
fraction of Poynting flux will. The fraction of a Poynting flux energy needed to significantly affect these results is yet unclear 
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